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Abstract. After presenting a brief summary of functional approaches to QCD at vanishing temperatures and 
densities the application of QCD Green's functions at non-vanishing temperature and vanishing density is dis- 
cussed. It is pointed out in which way the infrared behavior of the gluon propagator reflects the (de-)confinement 
transition. Numerical results for the quark propagator are given thereby verifying the relation between (de-)- 
confinement and dynamical chiral symmetry breaking (restoration). Last but not least some results of Dyson- 
Schwinger equations for the color-superconducting phase at large densities are shown. 



1 Introduction: Why Green's functions? 

The QCD Green's functions, and thereby especially the 
ones in the Landau gauge, have been in the focus of many 
recent investigations. In principle, they provide a complete 
description of the Strong Interaction. This implies that they 
also embody the non-perturbative phenomena of QCD, a- 
mongst them most prominently confinement, dynamical 
chiral symmetry breaking, and the axial anomaly. On the 
other hand, they serve as input into hadron phenomenology 
based on bound-state equations: mesons and their prop- 
erties are studied from solutions of Bethe-Salpeter, and 
baryons from the solutions of covariant Faddeev equations^ 

Of course, this raises immediately the question whether 
Green's functions are suitable to study the properties of 
the different QCD phases and the phase diagram of QCD. 
Although this question will likeley be answered affirma- 
tively, the investigations, as described in the following, are 
definitely at an early stage and require several improve- 
ments before they can be considered conclusive. Neverthe- 
less, they are certainly not only very promising but also the 
best access to map out the QCD phase diagram. 

Before continuing with a brief summary of functional 
approaches to QCD at vanishing temperatures and densi- 
ties, we want to point out that studies of the QCD phase 
diagram within functional methods have been discussed 
from a broader perspective than here in ^ (see also ref- 
erences therein). 



2 Functional Approaches to QCD 

As stated above, QCD Green's functions are a promising 
powerful tool. However, in order to determine them in con- 
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tinuum quantum field theory, in principle, an infinite hier- 
archy of coupled complicated equations has to be solved. 
As for QCD, we are of course most interested in their in- 
frared behaviour, i.e., in the Green's functions in the stron- 
gly-interacting domain. To this end it turns out that restrict- 
ing oneself to the primitively divergent Green's functions, 
which are the most interesting ones, provides an appro- 
priate starting point (see e.g. |4] and references therein). 
Even in the Landau gauge with the least number of prim- 
itively divergent functions a formidable task is left: there 
are seven primitively divergent functions while in Landau 
gauge Yang-Mills theory there are still five primitively di- 
vergent functions, namely the gluon and ghost propagators 
as well as the 3-gluon, 4-gluon and the gluon-ghost ver- 
tices. Including quarks, one has to consider in addition the 
quark propagator and the quark-gluon vertex. 

The functional approaches include methods which are 
based on the Dyson-Schwinger Equations (DSEs), the Func- 
tional Renormalisation Group Equations (ERG), and the 
«-Particle Irreducible Actions (nPI). Hereby, the DSEs, be- 
ing the equations of motion for the Green's functions, are 
derived straightforwardly from the generating functional. 
The ERG is based on the idea of employing energy-momen- 
tum cutoffs to "integrate out" the high-momentum modes 
and the nPI are effective actions allow to derive symmetry- 
preserving equations for the Green's functions. 

The advantages of functional methods are: 

- It is straightforward to implement chiral fermions and 
Goldstone's theorem. 

- There exist analytical infrared solutions. 

- Hadrons are described in terms of their fundamental 
substructure. 

- There is no sign problem at non-vanishing chemical 
potential. 

Compared to the most widely employed non-perturbative 
first-principle method, the lattice gauge theory, these ad- 
vantages are clear benefits. But these do not come for free, 
because, on the other hand, functional methods miss the 
main advantages of lattice gauge theory: 
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Fig. 1. The Landau gauge gluon renormalization function ob- 
tained from different calculations, see text for details. 

- There are no truncations. 

- Manifest gauge invariance. 

This makes plain that, whenever possible, an intelligent 
combinations of methods will provide the most reliable re- 
sults. 



3 What do we know for r = and // = 0? 

3.1 Infrared Structure of Landau gauge Yang-Mills 
theory 

In order to demonstrate the progress of functional meth- 
ods over the last decade it is interesting to compare dif- 
ferent steps in the calculation of the Landau gauge gluon 
propagator with corresponding lattice results. Many lattice 
calculations of the gluon propagator are available, and it 
would be beyond the scope of these proceedings to cite 
only the major steps in improving on these calculations. 
In Fig. [T] we have chosen the lattice results for the gluon 
renormalization function (see also Eq. ([TJ) of Ref. [5] for 
comparison, but any other recent calculation would equally 
well be suitable. The lowest lying line displays the results 
of Ref. fS] in which the coupled DSEs for the gluon and 
ghost propagators have been solved employing some ap- 
proximative treatment of the involved angular integrals. 
The second lowest curve is from Ref. IJI in which the an- 
gular integrals could be treated without any approximation. 
The curve above is the result of Ref. [8 1 where the ERG has 
been used to solve for the propagators. The highest lying 
curve which becomes also closest to the lattice data is a re- 
sult obtained by Jan Pawlowski based on the work of Ref. 
Il9l where DSE and ERG methods have been combinedQ 

How has this success of functional methods been a- 
chieved? To this end one should note that the ghost-gluon 
vertex becomes bare in the limit of vanishing ghost mo- 
mentum due to the transversality of the gluon propaga- 
tor This has been noticed already 40 years ago 111 Oil and 
has been tested in several non-perturbative investigations 

^ We thank Jan Pawlowski for discussing with us his work 
in preparation and Lorenz von Smekal for a compilation of 
the data (see also his talk under http://www.thphys.uni-heidel- 
berg. de/~ smp/view/Main/Delta 1 OProgramme) . 



IIIIII2II3I . With the knowledge that the ghost-gluon vertex 
stays regular in the infrared one can now analyze the ghost 
propagator DSE. Power-law like ansatze for the gluon and 
ghost renormalization functions lead to a relation between 
the infrared exponents ||6l (this kind of relations are nowa- 
days known as scaling relations). Denoting by Z{p~) and 
G{p^) the gluon and the ghost renormalization function, 
respectively, one obtains for — > 

z(/) ~ (p^f, G(p^) ~ {p^r, (1) 

in terms of the ghost propagator infrared exponent k. As it 
can be shown that 0.5 < k < 1 lfT4l the gluon propagator 
is infrared suppressed and the ghost propagator is infrared 
enhanced^ 

This type of infrared analysis of DSEs can be extended 
to all Yang-Mills vertex functions 1T7]T810 Employing in 
addition the ERG equations, and requiring that these two, 
seemingly different, towers of equations have to give iden- 
tical Green's functions, allows one to restrict the type of so- 
lutions strongly. There is one unique scaling solution with 
power laws for the Green's functions 1201211 and a one- 
parameter family of solutions, the so-called decoupling so- 
lutions. At this point some notes are in order Quite some 
time ago, in the Coulomb gauge, a similar situation has 
been found in variational approaches Il22ll : A family of 
infrared trivial solution possesses as an endpoint a criti- 
cal solution characterized by some infrared power laws. 
As it concerns the Landau gauge: Numerical solutions of 
the decoupling type (there called "massive solution") have 
been published in II23I24II and references therein. A recent 
detailed description and comparison of these two type of 
solutions has been given in Ref. 19i0Eirstly, one has to 
note that the scaling solution respects the BRST symmetry 
whereas the decoupling solutions break this symmetry |)9j. 
Secondly, almost all lattice calculations of the gluon prop- 
agator favor a decoupling solution. But lattice studies at 
strong coupling 12712812911 reveal the existence of a regime 
where the scaling relation between the gluon and the ghost 
propagator is fulfilled, and the corresponding infrared ex- 
ponent K is very close to the value determined in truncated 
continuum studies with k - 0.595. A potential resolution 
of this puzzle has been offered recently in Ref. 1301 : The 
infrared behaviour of the Green's function depends on the 
non-perturbative completion of the gauge. 

At this point it has to be emphasized that this discus- 
sion is of fundamental importance (especially with respect 
to the role of the BRST symmetry in a non-perturbative 



^ It is an interesting side remark that with this behaviour the 
Kugo-Ojima confinement criterion, the Oehme-Zimmermann 
superconvergence relation, and the Gribov-Zwanziger horizon 
condition are fulfilled, see, e.g., the review 1151 for a detailed 
discussion. 

Note that a recently pubhshed MATHEMATIC A code can be 
used to simplify significantly the derivation of the corresponding 
DSEs |T9). 

^ The infrared analysis for both type of solutions is described 
in Ref s. I25l26l . 



Hot and Cold Baryonic Matter - HCBM 2010 



approach) but the difference between the two types of so- 
lutions is phenomenologically irrelevant^ 

For completeness we will give the infrared behaviour 
of all one-particle irreducible Green's functions in the scal- 
ing solution in the simplified case with only one external 
scale p~ — > 0. For a function with n external ghost and 
antighost as well as m gluon legs one has: 

r"'"\p^) ~ {p^)^"-"'^' . (2) 

Note that this solution fulfills all DSEs, all FRG equations 
and all Slavnov-Taylor identities. In addition, it verifies the 
hypothesis of infrared ghost dominance [311 and leads to 
infrared diverging 3- and 4-gluon vertex functions. 

As a further side remark we want to add that the gluon 
propagator violates positivity 021331 which is a property 
related to the confinement of transverse gluons. Further- 
more, in Ref. Il32l an analytic structure for the gluon prop- 
agator has been proposed in which the gluon propagator is 
analytic in the complex p^-plane except the space-like real 
half-axis. 



3.2 Quarks: Confinement vs. DxSB and the Ua(1) 
anomaly 

Independent of the type of the solutions - scaling or decou- 
pling - the gluon propagator is infrared suppressed, and 
therefore quark confinement cannot be generated by any 
type of gluon exchange together with an infrared-bounded 
quark-gluon vertex. To add quarks one considers the DSE 
for the quark propagator together with the one for the quark- 
gluon vertex in a self-consistent way 041351 . Hereby, a 
significant difference of the quarks compared to the Yang- 
Mills fields has to be taken into account: Since they pos- 
sess a mass, and since dynamical chiral symmetry break- 
ing (D^'SB) does occur, the quark propagator will always 
approach a constant in the infrared. 

In general, the fully dressed quark-gluon vertex can 
be expanded in twelve linearly independent Dirac tensors. 
Half of the coefficient functions would vanish if chiral sym- 
metry were realized in the Wigner-Weyl mode. From a 
solution of the Dyson-Schwinger equations we infer that 
these "scalar" structures are, in the chiral limit, gener- 
ated non-perturbatively together with the dynamical quark 
mass function in a self-consistent fashion. Thus, dynami- 
cal chiral symmetry breaking manifests itself not only in 
the propagator but also in the quark-gluon vertex. 

Based on the Yang-Mills scaling solution one obtains 
from an infrared analysis an infrared divergent solution 
for the quark-gluon vertex such that the Dirac vector and 
"scalar" components of this vertex are infrared divergent 
with an exponent -k - \ if either all momenta Il34l or the 
gluon momentum vanish fi5\. A numerical solution of a 
truncated set of DSEs confirms this infrared behaviour. The 
essential components, to obtain this solution, are the scalar 
Dirac amplitudes of the quark-gluon vertex and the scalar 
part of the quark propagator. Both are only present when 

* Maybe with the exception of how to describe the axial 
anomaly with Green's functions, see below. 



chiral symmetry is broken, either explicitely or dynami- 
cally. 

This self-consistent quark propagator and quark-gluon 
vertex solution relates to quark confinement via the anoma- 
lous infrared exponent of the four-quark function. The static 
quark potential can be calculated from this four-quark one- 
particle irreducible Green function, which behaves like 
ip^Y' for — » due to the infrared enhancement of the 
quark-gluon vertex for vanishing gluon momentum. Using 
the well-known relation for a function F oc {p^Y^ one gets 

r 7^^'(/ = 0,p)e'P'- ~ |r| (3) 

for the static quark-antiquark potential V(r). Therefore the 
infrared divergence of the quark-gluon vertex, as found 
in the scaling solution of the coupled system of DSEs, 
the vertex overcompensates the infared suppression of the 
gluon propagator such that one obtains a linearly rising po- 
tential. 

Generally, if an approach is complete it has to describe 
all phenomena of a given theory. This especially implies 
that also the axial Ua(1) anomaly should be described with 
the Green's functions, or in other words, that the QCD 
Green's functions incorporate the topological susceptibil- 
ity of the QCD vacuum. To this end we would like to point 
out that the infrared divergence of the quark-gluon vertex, 
mentioned above, indeed generates an 77' mass and there- 
fore describes the axial anomaly |[36l . Following an old 
idea |371 which attributed the 77' mass to the infrared slav- 
ery of quarks it was shown that the infrared divergence of 
the quark four-point function in the soft limit is exactly 
the one which is needed to find a non-vanishing rj' mass. 
However, the appearence of the correct infrared singular- 
ity in single Feynman diagrams is the case only for the 
scaling solution. To our best knowledge, it is up to now 
unknown how the axial anomaly is encoded in the elemen- 
tary Green's functions of the decoupling solution as only 
a resummation of infinitely many diagrams will be able to 
describe the anomaly when employing these solutions for 
the Green's functions. 



4 Phases of strongly-interacting matter: 
How to go to r ^ and iui^0 7 

The exploration of the QCD phase diagram, in particular 
the higher baryon-density regime, within functional ap- 
proaches is a very active field of research and of great im- 
portance for a deeper understanding of the experimental 
data of the running and planned heavy-ion programs. As 
already argued, in the last years a fruitful interaction be- 
tween the different functional methods has lead to a largely 
quantitative understanding of QCD at vanishing temper- 
ature and density while an understanding of the confine- 
ment mechanism and its relation to the Dji'SB is not yet 
settled. At finite temperature and chemical potential the 
situation is even much less clear. On the one hand, with 
the help of functional techniques results for the pure Yang- 
Mills sector of QCD at finite temperature as well as for 
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the hadronic sector of QCD could be obtained but on the 
other hand a full QCD investigation is hampered by the fact 
that the gauge sector, i.e. the (de-)confinement transition, 
is not fully resolved. Furthermore, for full QCD with dy- 
namical quarks the (de-)confinement transition is expected 
to be a smooth crossover as the underlying center sym- 
metry of the SU(N(.) gauge group is explicitly broken by 
the quarks. Similarly, the nature of the chiral phase tran- 
sition depends, among other quantities, on the values of 
the current quark masses which again break chiral symme- 
try explicitly. It is an important observation that for small 
chemical potentials both phase transitions lie with respect 
to the temperature remarkably close to each other This is 
interesting since the (de-)confinement transition is driven 
by the gluodynamics while the chiral transition is govern 
by strongly-interacting quarks. The deeper understanding 
of this interrelation is a primary focus of the application of 
functional approaches, see also 031161 . 

A first step towards full QCD with functional methods 
at non-vanishing temperature and densities is done in the 
framework of eff'ective models which are generally con- 
structed from non-perturbative Yang-Mills effective poten- 
tials and hadronic potentials. Pure Yang-Mills theory cor- 
responds to the heavy-quark limit of QCD where the ex- 
pectation value of the Polyakov loop serves as an order pa- 
rameter for the (de-)confinement transition. In these mod- 
els some information about the confining glue sector of 
QCD is incorporated in an effective Polyakov-loop poten- 
tial that is extracted from pure Yang-Mills lattice simula- 
tions at vanishing chemical potential. This sector is effec- 
tively combined with the chiral quark(-meson) matter sec- 
tor leading to the Polyakov-NJL-type effective chiral mod- 
els such as the PNJL or PQM models, e.g., O8I39I40I . The 
matter sector of these models has been studied also beyond 
mean-field level by taking into account the quark-meson 
quantum fluctuations within a FRG approach, for recent 
studies see Refs 1411421431 . 

The most difficult problem within these model investi- 
gations is the inclusion of the quantum back-reaction of the 
matter sector to the gluonic sector. This problem has been 
attacked in Ref. |[39J which leads to a flavor and chem- 
ical potential dependence of the (de-)confinement transi- 
tion temperature. This first perturbative estimate of Ref. 
||39l has been confirmed by first-principle QCD calcula- 
tions within the FRG approach at real and imaginary chem- 
ical potentials 1161 and also by constraining PNJL model 
results with those in the statistical model, see 1441 for a 
recent review. This akeady demonstrates that even with 
QCD-based effective models combined with functional ap- 
proaches valuable information about the QCD phase di- 
agram can be provided and certain scenarios can be ex- 
cluded. As an example see 1431 concerning the expected 
focusing of isentropic trajectories in the vicinity of a criti- 
cal endpoint in the phase diagram. In this sense these mod- 
els can be understood as controlled approximations to full 
dynamical QCD. 

Nevertheless, treating the QCD Green's functions di- 
rectly is, of course, very desireable. In the next two subsec- 
tions two types of investigations aiming at a direct calcula- 



tion of the QCD Green's function within functional meth- 
ods will be described. 



4.1 T ?i and yu = case 

At r the heat bath provides a prefered rest frame. Ac- 
cordingly, the Green's functions will not only depend on 
the Lorentz invariant combinations of energy and momenta. 
This, by itself, makes the solution of functional equations 
much more involved. In addition, the number of functions 
which specify uniquely a given Green's function increases. 
An exception is the ghost propagator which still can be de- 
scribed by only one function. However, for the gluon prop- 
agator the one renormalization function spUts into two, 

^ou,o„ ^ M^pl^^p) ^ ^P^P)^ (4) 

according to transverse and longitudinal projections to the 
heat bath. For further discussions it is important to note 
that the corresponding contributions are essentially chro- 
momagnetic and chromoelectric in nature. For the quark 
propagator also the renormalization function splits into two 
such that the quark propagator is described by three func- 
tionfl 

1 

Sip) = . (5) 

-ijiPiAip) - iy4CjpC{p) + B(p) 

AtT-0 one has of course Zt = = Z and C = A 
which are functions of p^ only with p - (aip, p). At T 
we deal with six propagator functions depending on and 
Wp. At r — > oo one obtains a three-dimensional Yang-Mills 
theory plus a scalar field originating from the A4. 

At any T the chromomagnetic part of the gluon prop- 
agator keeps positivity violating I46I47I48 I . This clearly 
indicates a kind of partial gluon confinement at any tem- 
perature (a more detailed discussion may be found in [491 
and references therein). The results of these investigations 
relate also to the non-perturbative origin of chromomag- 
netic screening, whatever the temperature is. This is not at 
all surprising because three-dimensional Yang-Mills the- 
ory is confining (and thus the infrared behaviour of the 
gluon propagator is genuinely non-perturbative), one has 
an area law for the spatial Wilson loop, and the Coulomb 
string tension is non-vanishing at any T ll50l . 

To calculate the quark propagator we solve its DSE 
with the finite-temperature gluon propagator and quark- 
gluon vertex as input. In our corresponding first calcula- 
tions we follow Ref. 1511 where a fit to the lattice gluon 
propagator and an ansatz for the quark-gluon vertex has 
been used. As already stated the chromomagnetic part of 
the gluon propagator is not at all influenced by the phase 
transition. On the other hand, the chromoelectric part may 
serve as an order parameter. In Fig.|2]one sees very clearly 

' By symmetries a structure of the type ■yjpjy4a>pD(p) would 
be allowed which, however, has turned out so far in all investi- 
gated cases as negligibly small. 
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Fig. 2. Displayed is the Landau gauge chromoelectric screening 
mass scale with a fit to the critical exponent near the critical point. 
The data points are taken from Ref. 1511 . 
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Fig. 3. The Landau gauge scalar quark function B{p-) (see Eq. 
^) for dilferent temperatures in the chiral limit (solid lines) and 
for a typical light-quark current quark mass (dashed lines). 
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Fig. 4. The Landau gauge quark renoiTnalization functions A(p^) 
and C(p^) (see Eq. l|5]l) for different temperatures in the chiral 
limit. 



the drastice change in the infrared value of the chromo- 
electric gluon propgator, Dl{0), at T - Tc =277 MeV. The 
inverse of the square root of Dl{0) is the chromoelectric 
screening mass scale, and from the fit in Fig.|2]we obtain a 
critical exponent of approximately 0.53. 

The resulting quark propagator functions from the cor- 
responding DSE are displayed in Figs.[3]andj4] The scalar 
function B{p^) is the one reflecting D^^SBIj Correspond- 
ingly, in the chiral limit it vanishes identically above Tc 
thereby signaling the chiral phase transition. Of course, 
this transition becomes a crossover for non-vanishing cur- 
rent quark mass, and the function B is slowly varying and 
mostly driven by the current quark mass. A not anticipated 
behaviour is the fact that B{Q) raises with temperature up 
to Tc- But Fig.|2]makes evident why this happens: The in- 
frared value of the chromoelectric part of the gluon prop- 
agator increases with T up to before it then sharply de- 
creases. 

Putting as usual antiperiodic boundary conditions for 
quarks, the chiral-limit quark condensate (being essentially 
the functional trace over the quark propagator) serves as an 
order parameter However, playing with the boundary con- 
ditions allows to link confinement with spectral properties 
of the Dirac operator [|52l . CoiTespondingly, in Ref. | i3T| 
also the dual quark condensate and the dressed Polykaov 
loop II52I53I54I have been calculated^ The trick consists 
in setting t/(l)-valued boundary conditions for the quark 
field in "temporal" direction: 



q(x,/3^ l/T)^e"fq(x,0), 



(6) 



which amounts to an effective shift of Matsubara frequen- 
cies, 

— 2nT{n + ipjl-Ti). (7) 

The generalized condensate (^q)^, at fixed is the corre- 
sponding expectation value of the Dirac operator, and it is 
expandable in a geometric series containing loops of link 
variables with increasing winding number. One can then 
define the dual (Gattringer) condensate as 



Jo 27r 



(8) 



where v — \ projects out winding-number-one loops, the 
so-called dressed Polyakov loops Il52l . It is an order pa- 
rameter for center symmetry breaking and confinement. 

The aim of our future investigations is now to calcu- 
late the dual quark condensate and the dressed Polyakov 
loop with an improved input. Substituting the parameter- 
ized gluon propagator by recent results of functional equa- 
tions is hereby more a matter of convenience and numer- 
ical precision than of fundamental progress. However, an 
(at least partially) self-consistent determination of the quark 
propagator and the quark-gluon vertex together, generaliz- 
ing the work of Ref. 1351 to non-vanishing temperatures. 



* For conceptional clarity we prefer to discuss the function 
S(p^) and not one of the mass functions M(p^) = B(p^)/A{p^) 
or M(p^) = B(p^)/C(p^), respectively. 

' For corresponding T = lattice and functional equations see 
Refs. Il55l56l . 
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would be a decisive step towards a first-principle calcu- 
lation of finite-temperature quark properties in continuum 
QCD. 

At the end of this subsection we want to point out that 
one can define a Polyakov-loop related confinement cri- 
terium from infrared exponents Il57ll . Together with an FRG 
calculation of the Polyakov potential this allows one to cal- 
culate the transition temperature. In pure Yang-Mills the- 
ory one finds hereby for SU(2) a second-order phase tran- 
sition lying in the Ising universality class and a first-order 
transition for all SU(Nc >3), Sp(2) and for the E(7) group 



4.2 Color superconducting phase 

Of course, it is quite obvious to exploit the absence of the 
sign problem in functional approaches to extend the stud- 
ies to non-vanishing density. Nevertheless, the introduc- 
tion of a chemical potential into functional equations leads 
to a significant complication of the quark propagator pa- 
rameterization, usually done in the Nambu-Gorkov formal- 
ism then, see e.g. Refs. 1591601611 and references therein. 
It has to be noted that as long as one uses an Abelian- 
type model for the quark-gluon vertex these difficulties 
are merely technical obstacles which can be overcome by 
combining computer algebra with numerical techniques. 
Nevertheless, including medium modifications by quarks 
(as done e.g. in Refs. 1591601611 ) and employing the re- 
sults of functional equations or a combined phenomeno- 
logical and lattice fit to the gluon propagator (as e.g. in 
Ref. 1621 ) one obtains a quite astonishing but robust result: 
Restricting oneself to translationally invariant phases one 
concludes that for a realistic renormalized strange quark 
mass for all chemical potentials above the phase transition, 
quark matter is in the colour-flavour locked colour-super- 
conducting phase. In Fig. |5]a resulting phase diagram in 
the plane of renormalized strange quark mass and quark 
chemical potential is given^The shaded band in the figure 
indicates the value of the strange quark mass as given by 
the Particle Data Group, the different lines are phase sep- 
aration lines from colour-flavour-locked to gapless colour- 
flavour-locked (lowest line), from gapless colour-flavour- 
locked to the uSC phase (second-lowest line), and so on. 

In addition, these calculations demonstrate that there 
are huge deviations of the gap functions as compared to 
those extrapolated from the weak-coupling result up to che- 
mical potentials of the order of several GeV. Even at such 
large chemical potentials perturbation theory quantitatively 
fails. Furthermore, the light quarks screen interaction also 
in the strange quark sector, an effect, which is not seen in 
most models. 

A calculation of the dressed Polyakov loop at finite 
chemical potential will certainly shed more light on the 
nature of the finite-density phase transition, and therefore a 
corresponding generalization of the work described in Ref. 
II59I60I6III is in progress Il63l . However, it has to be stated 

To obtain this figure electric and colour neutrality have been 
imposed, but this is, however, not decisive for the main result. 
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Fig. 5. Color superconducting phases in the - /i plane. 

clearly here that the main aspects of the QCD phase tran- 
sition at finite density will not be uncovered by functional 
methods if we will not get more information on the cou- 
pling of quarks to gluons also in this region of the phase di- 
agram. In the case of a DSE study this especially includes 
more knowledge about the quark-gluon vertex 

5 Summary and Outlook: What may we 
expect? 

The last years have seen a quite dramatic increase in our 
knowledge on the Landau gauge QCD Green's functions 
at vanishing temperatures and densities. It has become evi- 
dent that the propagation of transverse gluons is infrared 
suppressed, and that positivity is violated for transverse 
gluons. It has turned out that ghosts are the long-range "de- 
grees of freedom", and concentrating on the main point one 
may say that in Landau gauge QCD, ghosts and the restric- 
tion to the Gribov horizon {i.e. the correlations introduced 
by gauge-fixing uniquely) are related to the origin of con- 
finement (|3lll5l . 

It cannot be emphasized enough that dynamical chiral 
symmetry breaking does not only lead to the generation 
of constituent quark masses but also to scalar-type cou- 
plings between quarks and gluons. There are indications 
that the related scalar confinement potential is even larger 
than the vector component ll35ll . We also remark that the 
axial anomaly is encrypted in the infrared behaviour of the 
Green's functions in the scaling solution 1361 . 

Having set the stage, one is ready to continue to ei- 
ther non-vanishing temperatures or/and densities. Firstly, 
one should note that the chromomagnetic part of the gluon 
propagator does not see the phase transition at all. In par- 
ticular, one has positivity violation (and thus partial gluon 
confinement) at any temperature. Secondly, the infrared re- 
gion of the chromoelectric part of the gluon propagator dis- 
plays very nicely the phase transition ISTl . Needless to say, 
that in the chiral limit the quark propagator changes also 
accordingly due to the chiral symmetry restoration at the 
critical temperature. 
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The results of the quark propagator DSE at non-vani- 
shing densities provide evidence that for a realistic renor- 
malized strange quark mass for all values of the chemical 
potentials above the phase transition, quark matter is in the 
colour-flavour locked colour-superconducting phase IISTI . 
However, here we have to keep in mind that these investi- 
gations are restricted to translationally invariant phases. 

As we know from recent studies within functional ap- 
proaches at r = and - Q the dressing of the quark- 
gluon vertex is of utter importance. Therefore we will con- 
centrate in future on investigations of this function at non- 
vanishing temperatures and densities. Hereby, the dual qua- 
rk condensate and the dressed Polyakov loop 155! will be 
important quanitities in the studies of the QCD phase di- 
agram. Although it may sound very ambitious today we 
want to emphasize that the goal is to calculate thermody- 
namic observables of QCD at all physically relevant tem- 
peratures and chemical potentials. 
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